Dynamical percolation transition in the Ising model studied using a pulsed magnetic 

field 
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We study the dynamical percolation transition of the geometrical clusters in the 2d Ising model 
when it is subjected to a pulsed field below the critical temperature. The critical exponents are 
independent of the temperature and pulse width and are different from the (static) percolation 
transition associated with the thermal transition. For a different model that belongs to the Ising 
universality class, the exponents are found to be same, confirming that the behavior is a common 
feature of the Ising class. These observations, along with a universal critical Binder cumulant value, 
characterises the dynamical percolation of the Ising universality class. 



I. INTRODUCTION 

The thermal phase transition of the Ising model have 
been extensively studied in the geometrical terms of (cor- 
related) percolation in the past. It was seen that the "ge- 
ometrical" clusters, consisting of nearest neighbor par- 
allel spins, undergo a percolation transition [ljtlj at the 
same critical temperature as the magnetisation in two di- 
mensions 0, @. In higher dimensions, however, the tran- 
sition point differs [l|, 0] • Even for the two dimensional 
case, although the critical points coincide, the critical be- 
havior does not @, Q. Later, by redefining the clusters 
suitably Q (the so called "physical" clusters), the critical 
point was made to coincide in all dimensions along with 
the critical behavior fiol ]. 

Although the temperature induced percolation transi- 
tion is well studied, very little is known about the per- 
colation transition induced by external magnetic field. 
It is well known [TTl - fl7j that pure Ising systems, below 
its static critical temperature (T®) can undergo a field- 
pulse induced magnetisation reversal transition. The 
field pulse is applied to the Ising system in the direc- 
tion opposite to that of the existing order for a finite 
time (At). Depending on the value of the field amplitude 
(h p ), temperature (T) and pulse duration (At), the sys- 
tem can eventually undergo a transition from one equilib- 
rium state (with magnetisation —mo) to the other (with 
magnetisation mo) after the field is withdrawn. Apart 
from its theoretical interests, it has experimental reali- 
sations [l8| as well as industrial applications, such as in 
recording and switching industry (for review see 

As there is a percolation transition associated with 
the thermal one, here also we expect a similar perco- 
lation transition. In the following study, we investigate 
this percolation transition of the "geometrical" clusters 



at different temperatures, field pulse amplitude and field 
pulse widths. We find the critical exponents associated 
with the percolation transition to be significantly differ- 
ent from the percolation transition of the "geometrical" 
clusters in the thermal (or static) transition. These expo- 
nents remain invariant with change in temperature and 
field pulse width. We have also measured the critical 
Binder cumulant value 
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of the percolation transition 
order parameter. This quantity is not only useful in lo- 
cating the transition point but also its value at the critical 
point (which is independent of system size) indicates the 
universality class to which the transition belongs. It is 
seen that this value is different from the thermal coun- 
terpart and is very robust in the sense that it remains 
unchanged with changes in temperature and field pulse 
width. 

It was shown before [22|, [23| for the static transition 
that the percolation exponents of the geometrical clus- 
ters have same values for the models belonging to the 
same universality class (the case of Ising and Z(3) sym- 
metric models were studied). For all models belonging to 
Ising class, therefore, universal exponents were found. In 
our study of dynamical percolation transition also, if the 
exponents have to characterise the dynamical percolation 
transition of the Ising universality class, their values have 
to be same for all models belonging to this class. We have 
checked for a different model (in Ising class) that the ex- 
ponents and the critical Binder cumulant value are the 
same. 

For the magnetisation reversal transition also, it was 
claimed [ill, Eij that the universality class was different 
from the static one. In that case, although the exponents 
were robust with respect to changes in temperature and 
pulse width, we find that the critical Binder cumulant 
value changes continuously with temperature, leaving the 
claim incomplete. However, as the associated percolation 
transition (occurring at the same critical point within our 
numerical accuracy) clearly belongs to a different univer- 
sality class, one can probably say that the claim of dif- 
ferent universality class there was actually true. And the 
non-universal critical Binder cumulant value can be put 
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as yet another example of some special cases, which is 
known to occur also for changes in boundary condition, 
lattice shape, anisotropic effect etc. [24T - l32j . 

The rest of the paper is organised as follows, in sec. II 
we define the model and the percolation transition of the 
"geometrical" clusters. In sec. Ill we measure the critical 
exponent and Binder cumulant values for different tem- 
peratures and pulse widths and compare them with the 
magnetisation reversal transition. In sec. IV wc discuss 
our main results and conclude. 



II. FIELD PULSE INDUCED PERCOLATION 
TRANSITION 

A. Pure Ising model 

Here we have studied the percolation transition induced 
by an external magnetic field pulse in the case of pure 
two dimensional Ising model with nearest neighbour in- 
teraction. The Hamiltonian of the system reads 

H = -J^J ij S i S j -h(t)J2S i , (1) 

{ij} » 

where Si's (i = 1,2, ... , N) are Ising spins on lattice site i 
and Jij is the cooperative interaction between the near- 
est neighbor spins. Since we have considered a square 
lattice, the static critical point T c ° = 2.269... (in units of 
J/fcs). Periodic boundary condition is used along both 
directions. We keep the temperature T of the system 
below this critical value such that there is a spontaneous 
equilibrium magnetisation. The time dependent external 
magnetic field h(t) is applied for a finite duration in the 
direction opposite to that of the equilibrium magnetisa- 
tion. Although spatially uniform, the time dependence 
of this applied field is as follows: 

Ht) = l V *o<i<*o + At 
v ' 10, otherwise. v ' 

In a given Monte Carlo (MC) update a spin is randomly 
selected and the energy difference (AE) caused for flip- 
ping it is calculated. The flip is made if a random number 
is less than exp(— AE/T). L 2 such updates make a sin- 
gle MC step. Here wc take the geometrical definition of 
the clusters where the clusters are formed from parallel 
spins and each spin in a cluster must have at least one of 
its neighbors parallel to it. At t = to, there is a sponta- 
neous magnetisation (say, in the "downward" direction) 
in the system and one spanning cluster (of "down" spins) 
is present. Then magnetic field is applied in the opposite 
direction ("upward") for a finite duration (At) and the 
absolute value of the magnetisation as well as the per- 
colation order parameter P max = Sl/L 2 (where Sl is 
the size of the largest cluster and L is the linear size of 
the system) decreases with time in that duration. For a 
particular combination of T and At at some hp(T,At), 
the system undergoes a percolation transition as more 



and more spins flip upward and the previously percolat- 
ing cluster (of down spins) becomes non-percolating (sec 
Fig. [T|). In the following, we have studied the critical 
behavior of this percolation transition at different points 
in the h p — At phase boundary (details are shown for 
one point only) which coincides with (within our numer- 
ical accuracy) the phase boundary obtained before [ll[ 
in the magnetisation reversal transition (where the order 
parameter was defined [Tlj as the magnetisation m w at 
which the field is withdrawn). 
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FIG. 1: (Color online) Variations of the order 
parameters in the two transitions (a) withdrawal point 
magnetisation (m w ) and (b) the probability of largest 
cluster (Pmax) with the field pulse amplitude (h p ) in 
the case of pure Ising model, where L = 100, T = 1.0 
and At = 4. 



B. Ising model with diagonal second neighbor 
frustration 

To verify if the transition mentioned above characterise 
the dynamical percolation of the Ising universality class, 
here we study the dynamical percolation in a different 
model that belongs to the Ising universality class. This 
model of Ising spins has nearest neighbor (NN) ferro- 
magnetic and (diagonal) next nearest neighbor (NNN) 
(weaker) antiferromagnetic coupling 

H = -J x ]T SiSj + S i S i + M*) Si - ( 3 ) 

NN NNN i 

We study this model for J2/J1 = 1/10 (as was chosen 
in [22|]) and h(t) following Eqn. (2). The static (when 
h(t) = 0) critical temperature for the order-disorder tran- 
sition in this model is T c 1.945 [22j in the same unit as 
before. We perform the following study below this critical 
temperature (at T = 1.0). At that temperature, the sys- 
tem is in ferromagnetic order. As before, a magnetic field 



3 



pulse is applied in the direction opposite to this order for 
a finite duration. For a fixed value of temperature and 
pulse width (4 MC time steps), by increasing the field 
pulse amplitude h pi the previously spanning cluster of 
parallel spins becomes non-spanning and also the system 
undergoes a magnetisation reversal transition (Fig 2). 
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FIG. 2: (Color online) Variations of the order 
parameters of magnetisation reversal (m w ) and the 
percolation transition (P ma x) with the field pulse 
amplitude (h p ), where L = 100, T = 1.0 and At = 4 in 
the case of Ising model with next nearest neighbor 
frustration. 

It should be mentioned here that when there is no 
external magnetic field, the thermal transition and the 
percolation transition of the geometrical clusters can be 
shown 0, Q to occur at the same critical point in two di- 
mensions. In our case, however, the transition is driven 
by field which is time dependent. Hence, the evidence of 
coincidence in this case is entirely numerical. 



III. CRITICAL BEHAVIOR: FINITE SIZE 
SCALING 

The percolation transition is characterised by power-law 
variation of different quantities. The order parameter 
i.e., the relative size (P m ax) of the largest cluster varies 
as 



(4) 



The correlation length diverges near the percolation tran- 
sition point as 



£ ~ (K - M 



(5) 



where, hp is the critical field amplitude. 
The values of the critical exponents /3 and v specify the 
universality class of the transition. The other exponents 
can be obtained from scaling relations [H|. 
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FIG. 3: (Color online) Inset: P max ^ lv is plotted 
against the field amplitude (h p ) where T = 1.0 and 
At = 4 for pure Ising model. The curves for different 
system sizes (L = 100, 200, 400, 500) cross at 
h c p = 2.105 ± 0.005 giving the critical point and 
P/v = 0.20 ± 0.05. Front: The value obtained for /3/v is 
used in the y-axis and by tuning the value of 1/v all the 
curves for different system sizes were made to collapse 
on a single curve, thereby estimating 1/v = 0.85 ± 0.05. 
Error bars arc smaller than the symbol size. 



However, the exponents are not determined from these 
definitions due to finite size effects. The critical expo- 
nents are determined from the finite size scaling relations. 
For example, the order parameter is expected to follow 
the scaling form 



Pmax — L ^ J- L / {hp 



(G) 



where J 7 is a suitable scaling function. If we plot 
Pmax^ against h p for different system sizes, then by 
tuning P/v, all the curves can be made to cross at a 
single point. The field amplitude where this happens 
must be the critical field amplitude (hp). To estimate v, 
PmaxL®'" is to be plotted against (hp — h p )L 1 ' L ' and by 
tuning the curves are made to collapse, giving an 
accurate estimate of the exponent v. 

To further verify the critical point and the universality 
class, we study the reduced fourth order Binder cumulant 
of the order parameter, defined as (2ll | 



U = 1 - 



(P 4 ) 

\ max/ 

Q/p2 \: 
"\ max/ 



(7) 



where P ma x is the percolation order parameter (as de- 
fined before) and the angular brackets denote ensemble 
average. U — > | deep inside the ordered phase and U — >• 
in the disordered phase when the fluctuation is Gaussian. 

The Binder cumulant can also be used to find the cor- 
relation length exponent as it follows the scaling form 



U^U^p-hp)^), 
where hi is a suitable scaling function. 



(8) 
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A. Pure Ising model 

In Fig. [3] (inset) we plot P max L^' v against h p for a given 
temperature (T = 1.0) and field pulse width (At = 4), 
where the unit of time is measured by Monte Carlo Steps 
(one MCS being L 2 spin updates). The curves for dif- 
ferent system sizes cross at a given point when /3/v = 
0.20 ± 0.05 and the crossing point (h c p = 2.105 ± 0.005) 
gives the critical field amplitude. 

To estimate the exponent v, we plot for the same set, 
P m .axL p/u against (h c p - h v )L 1 l v (see Fig. EJ. Now by 
tuning the 1/v term in the X-axis only, all the plots 
are made to collapse on a single curve. This gives us the 
estimate of 1/v = 0.85±0.05. Again, the error bar is esti- 
mated by changing the tuning parameter upto the point 
when the scaling is visibly worsened. These estimates 
differ significantly from the corresponding static transi- 
tion values reported in Ref. [H] (fi s /v s = 0.052 ± 0.002, 
v s = 1.004 ±0.009). 
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also agrees well with our estimate of f3/v from finite size 
scaling (0.18 and 0.20 respectively). 




FIG. 4: (Color online) Size of the largest cluster for 
various system sizes 

(L = 100, 200, 300, 400, 500, 600, 700, 800, 900, 1000) at 
the critical field (h c p = 2.105) where T = 1.0 and At = 4 
for pure Ising model. The slope of the log-log plot gives 
the fractal dimension to be D = 1.82 ± 0.01. The other 
two curves show the variation of the same, below 
(hp = 2.07) and above (h p = 2.14) the critical field. 

It is known [33[ that the largest cluster at critical point 
has a fractal nature. One can find its fractal dimension 
from the relation Sl ~ L D , where Sl is the size of the 
largest cluster and D is its fractal dimension. The fractal 
dimension is related to the spatial dimension of the lattice 
through the critical exponents in the following way D = 
d — /3/v, where d is the spatial dimension. While for 
static transition in Ising model the fractal dimension of 
the critical droplet is 1.875 (for geometric cluster in static 
transition it is 1.947±0.002 jH), from the best fit (Fig.gJ 
of the largest cluster size with linear size in a log-log 
plot gives the value to be 1.82 ± 0.02 in our case, which 




FIG. 5: (Color online) Fourth order reduced Binder 
cumulant of the percolation order parameter (P max ) for 
different system sizes (L = 100, 200, 400, 500) are 
plotted for pure Ising model; crossing point determines 
the critical point (h p = 2.10 ± 0.01). The critical Binder 
cumulant value (U* = 0.52 ± 0.01) specifies the 
universality class. Inset shows the data collapse for 
same 1/v as obtained before (0.85). 

To check if the exponents vary with temperature 
or pulse width, we have scanned the parameter space 
(At, T) in the range T = 0.5 -)• 2.0 and At = 4 ->• 10 
with T = 0.5, 1.5, 2 and At = 4, 6, 10 (all combinations) 
apart from the combination (1.0,4) the details of which 
are shown above. We find no significant variation in the 
estimates of the critical exponents. 

The crossing point of the different curves for different 
system sizes gives the critical point (2.102±0.002), which 
is in good agreement with the previous estimation from 
finite size scaling (2.105 ± 0.005). The value of U at the 
critical point in this case is 0.52 ± 0.01 (see Fig. [5]). This 
value remains unchanged for different sets of (mentioned 
above) temperatures and field pulse width. In the inset 
of Fig. [5] the data collapse is shown. The same value of 
the exponent v is used, as was obtained before. 

Therefore, with the increase of temperature and/or 
field pulse width we find no crossover behavior in the 
percolation transition as was reported in 
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magnetisation transition (occurring within the range of 
the parameter space scanned here). There is indeed a 
single regime as far as the critical exponents and Binder 
cumulant values of the percolation transition are consid- 
ered. 



B. Ising model with second neighbor frustration 

The critical exponents of the dynamical percolation tran- 
sition in this case also are found using the finite size scal- 
ing analysis as is shown in the pure Ising case. The expo- 
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(hp C -V Ll/V 

FIG. 6: (Color online) Inset: P max L p/u is plotted 
against field amplitude (h p ) for the Ising model with 
frustration, keeping T = 1.0 and At = 4. The curves for 
different system sizes (L = 100, 200, 300, 400) cross at 
hp = 1.65 ± 0.01 giving the critical point for 
f}fv = 0.20 ± 0.05. Front: The value obtained for fijv is 
used in the y-axis and by tuning the value of 1/v all the 
curves for different system sizes were made to collapse 
on a single curve, thereby estimating l/v = 0.85 ± 0.05. 
Error bars arc smaller than the system size. 
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FIG. 7: (Color online) Fourth order reduced Binder 
cumulant of the percolation order parameter (P max ) for 
different system sizes (L = 100, 200, 300, 400) is plotted 
for the Ising model with frustration; crossing point 
determines the critical point = 1.65 ± 0.01). The 
critical Binder cumulant value (U* = 0.52 ± 0.01) 
specifies the universality class. Inset shows the data 
collapse for same l/v as obtained before (0.85). 



nents /3/V = 0.20 ± 0.05 and 1 jv = 0.85 ± 0.05 are found 
from data collapse shown in Fig. ©. 

The values compare very well with the estimates from 
the pure Ising model. Fig. (Q shows the behavior of 
the Binder cumulant of the percolation order parameter 
near the critical point and also its finite size scaling using 



cqn ©. The critical Binder Cumulant value (0.52±0.01) 
shows the expected agreement with the pure Ising case, as 
well as the previously estimated value of l/v is confirmed 
from the data collapse shown in the inset. 



IV. DISCUSSION AND CONCLUSION 

It is well studied that for pure Ising model in two dimen- 
sion, a phase transition in magnetisation and a percola- 
tion transition of the geometric clusters of nearest neigh- 
bor parallel spins occur at the same critical point. There 
have been successful attempts [13, [H| to show that this 
percolation transition is universal in the sense that its 
exponents are same for models in the same universality 
class, providing a new set of exponents that are charac- 
teristic features of Ising and other universality classes in 
terms of percolation. 

It is well known that pure Ising model can undergo 
a dynamical phase transition (magnetisation reversal) 
when subjected to competing pulsed magnetic field. This 
dynamical phase transition was claimed to belong to 
a different universality class than the static transition. 
However, we find that the critical Binder cumulant value 
reported in Ref [ll[ in not universal for that transition. 
It changes if the transition is studied, for example, in a 
different temperature or pulse duration etc. 

Here we have studied the same transition in terms of 
the percolation properties of the geometric clusters. The 
dynamical transition exponents were found using finite 
size scaling analysis. The values of these exponents re- 
ported here {fi/v = 0.20 ± 0.05 and u=1.2± 0.1), dif- 
fer significantly from those (I3 s /v s = 0.052 ± 0.002 and 
v s = 1.004 ± 0.009) obtained in Ref. [H] for the case of 
the static transition. Therefore, the dynamical percola- 
tion transition is in a different universality class than the 
static one. 

Further, the critical Binder cumulant value which was 
found to be non-universal for magnetisation (crossover 
behavior was reported), is found to be universal (doesn't 
changes with temperature and/or field pulse width, and 
microscopic details of the models in the same universal- 
ity class) when studied in terms of the percolation order 
parameter (within a reasonable range of the parameter 
space, see sec IIIB), which is the expected behavior (hav- 
ing unique value for a given universality class) of the crit- 
ical Binder cumulant. 

Also, to see if the exponents reported here indeed char- 
acterise the dynamical transition in Ising universality 
class, we have repeated the same analysis for another 
model (diagonal next nearest neighbor frustration) be- 
longing to the Ising universality class. It is found that 
the exponents obtained in this model matches very well 
with those of the pure Ising model. The critical Binder 
cumulant value also takes the expected universal value. 
Therefore it is concluded that the dynamical percolation 
behavior reported here is also a characteristic feature of 
any model belonging to the Ising universality class. 
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Finally, wc would like to point out that the values 
of the critical exponents obtained in our study has re- 
markable agreement with those obtained in [34[ in which 
the authors study percolation properties of randomly dis- 
tributed growing clusters using finite size scaling analysis. 
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